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The properties of uniformly rotating white dwarfs (RWDs) are analyzed within the frame¬ 
work of Newton’s gravity and general relativity. In both cases Hartle’s formalism is ap¬ 
plied to construct the internal and external solutions to the field equations. The white 
dwarf (WD) matter is described by the Chandrasekhar equation of state. The region of 
stability of RWDs is constructed taking into account the mass-shedding limit, inverse /3- 
decay instability, and the boundary established by the turning points of constant angular 
momentum J sequences which separates stable from secularly unstable configurations. 
We found the minimum rotation period ~ 0.28 s in both cases and maximum rotating 
masses ~ 1.534M0 and 1.516M0 for the Newtonian and general relativistic WDs, 
respectively. By using the turning point method we show that general relativistic WDs 
can indeed be axisymmetrically unstable whereas the Newtonian WDs are stable. 

Keywords-. Newtonian and general relativistic white dwarfs; maximum mass; minimum 
period; stability. 


1. Introduction 

Recently, equilibrium configurations of non-rotating (static) ^He, and ^®Fe 

white dwarfs (WDs) within general relativity (GR) have been constructed in Ref. 1. 
The white dwarf matter has been there described by the relativistic generalization 
of the Feynman-Metropolis-Teller (RFMT) equation of state (EOS) obtained by 
Rotondo et al.^ A new mass-radius relation that generalizes both the works of 
Chandrasekhar^ and Hamada & Salpeter"' has been there obtained, leading to a 
smaller maximum mass and a larger minimum radius with respect to the previous 
calculations. In addition, it has been shown how both GR and inverse /3-decay are 
relevant for the determination of the maximum stable mass of non-rotating WDs. 

It is therefore of interest to generalize the above results to the case of rotation. 
As a first attempt, we constructed in Ref. 5 general relativistic uniformly rotating 
WDs in the simplified case when microscopic Coulomb screening is neglected in the 
EOS, following the Chandrasekhar^ approximation by describing the matter as a 
locally uniform fluid of electrons and nuclei. The average molecular weight in the 
Chandrasekhar EOS is fi — AjZ^ where A is the mass number and Z is the number 
of protons in a nucleus. 

As a second attempt, in Ref. 6 we calculated the maximum mass of rotating 
^He, ^^C, ^®0 and ®®Ee WDs using the Salpeter^^ and the RFMT EOS. As a result 
we obtained there different maximum mass for different chemical composition of 
WD matter. 

As a third attempt, in Ref. 7 we investigated the stability of general relativistic 
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uniformly rotating "^He, and ^®Fe white dwarfs against secular and dynam¬ 

ical instabilities. We determine the maximum mass and minimum rotation period of 
stable white dwarfs depending on chemical composition of the white dwarf matter 
taking into account the Coulomb interactions as well as the nuclear interactions and 
the electroweak equilibrium at high densities, within the Salpeter^^ and relativistic 
Feynman-Metropolis-Teller EOS.^ 

We summarized all the above results in Ref. 8, where we analyzed in detail 
the stability of RWDs from both the microscopic and macroscopic point of view. 
Besides the inverse /3-decay instability, we also study the limits to the matter den¬ 
sity imposed by zero-temperature pycnonuclear fusion reactions using up-to-date 
theoretical models. 

In this work we show the importance of general relativity to study the stability 
of RWDs. Moreover one can restrict to the Newtonian case when estimating the 
minimum period of uniformly RWDs. In case of the maximum rotating mass and in 
the analyses of stability, general relativity indeed becomes relevant. We make use 
of Hartle’s formalism^^’^^ for the construction of rotating stars both in Newton’s 
gravity and GR. As a WD matter we consider the Chandrasekhar EOS with average 
molecular weight, /i = 2. 

The main parameters characterizing RWDs such as the total mass M, angular 
momentum J, quadrupole moment Q, equatorial radius Req and eccentricity e are 
obtained for a given EOS, from the matching procedure between the internal and 
external solutions of the field equations and the equations of hydrostatic equilib¬ 
rium. ^^^3 ]\jote, that the total mass is defined by M = + SM , where 

M-^=° = is the mass of a static white dwarf with the same central density as 
and 6M is the contribution to the mass due to rotation. 

The main objective of this work is to compare and contrast WDs both rotating 
and non-rotating in Newton’s gravity and general relativity. 


2. Stability criteria 

There are two main stability criteria for non-rotating WDs: the well-known gen¬ 
eral relativity (GR) instability and inverse /3-decay instability. The inverse /3-decay 
instability is related to the fact that with the increasing mass of WDs the central 
density will increase, hence the electrons will become ultra relativistic and their 
energy will increase until it reaches a certain value of the Fermi energy (depending 
on chemical composition of a WD matter) which allows the electrons to be captured 
by the protons inside the nuclei, forming neutrons. This process is called inverse /3- 
decay or neutronization and is written in a compact form p+e~ —> n + Ve, where p is 
the proton, e~ is the electron, n is the neutron and Ve is the electron neutrino. Such 
a process conserves the number of baryons creating neutron rich nuclei through e~ 
capture. Due to the fact that most of the pressure of the WD is supported by the 
electrons, the inverse /3-decay process leads to the pressure decrease of the WD by 
increasing the density which causes the instability of the star. 
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For rotating WDs there exist additional stability criteria. One of them is the 
mass shedding limit. When a star spins so rapidly the centrifugal forces will prevail 
over gravitational force and the star will start losing its matter on the equatorial 
plane.The procedure of estimating the mass shedding angular velocity in general 
relativity was developed in Friedman et. al.^® 

In addition to the mass shedding limit there are two more stability criteria called 
secular (axisymmetric) and dynamical instabilities. Usually in the literature the 
conditions T/W = 0.14, e = 0.81267 and T/W = 0.25, e = 0.95288 are adopted as 
the upper limit for the onset of the secular and dynamical instabilities,^"^ where TjW 
is the rotational/gravitational binding energy and e is the eccentricity of rotating 
objects. These values are absolute upper limits for the ideal case of Maclaurin 
spheroids^^ and in practice, they cannot be applied to objects like uniformly rotating 
WDs since they have smaller TjW and e with respect to the Maclaurin spheroids.^ 
However Friedman et. al.^^ developed the turning point method of axisymmetric 
(secular) instability for uniformly rotating stars. The main idea of the turning point 
method is that along a sequence of constant angular momentum, the condition 
{dM{p,J)/dp)j = 0 separates secularly stable and unstable stars. For J = 0 we 
have the condition for the maximum mass of non-rotating WDs dM{p)/dp — 0. 
Moreover it was shown in Friedman et. al.^^ that a sequence of uniformly rotating 
stars becomes dynamically unstable after it has become secularly unstable. That is, 
for fixed angular momentum the point of secular instability occurs at lower central 
density with respect to dynamical instability. 



M 

Wq 



— Static configuration in GR 

— Keplerian sequence in GR 

— Inverse ^S-decay inst. 

— Secular inst. in GR 

— Newtonian static configuration 

— Newtonian Keplerian sequence 

To L5 2^0 

l?eq[10^km] 


Fig. 1. Mass in solar mass units versus central density (left panel) and mass in solar mass units 
versus equatorial radius (right panel) of Newtonian and general relativistic WDs for the Chan¬ 
drasekhar EOS with /i = 2. Both the non-rotating case and the Keplerian sequence are shown. We 
have stopped the density, just for the sake of comparison, at the critical density for the onset of 
inverse /3-decay of “^He = 1.39x10^^ g/cm^. 


3. Results 

We found that along the Keplerian (the mass shedding) sequence there are no 
maxima of J =constant sequences. This means that along the Keplerian sequence 
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all uniformly rotating WDs are secularly, hence dynamically stable. For this reason 
we increase the central density all the way up to inverse ,5-decay density in order 
to define the minimum rotation period of WDs, since this is the only limitation 
which guarantees the stability before the trigger of gravitational collapse. Thus, the 
minimum period of WDs is determined along the Keplerian sequence at the critical 
central density for the inverse /3-decay instability. 




Fig. 2. Left panel: eccentricity e versus central density pc. Right panel: normalized quadrupole 
moment versus versus central density pc, where and R are the static mass and 

static radius, respectively. The Keplerian sequences are considered only. The solid red line is in 
Newton’s gravity and the solid black line is in general relativity. 


1.5 

. ::: 

1.5 


1.4 

M 


1.4 

M 


TT 1-3 

//' 

TT 1.3 

' '' 

Mq 

// 

Mq 


1.2 

// — Static configuration of this work 

\ 9 

— Static configuration of this work ^ \ \ 

// — Keplerian sequence of this work 


— Keplerian sequence of this work ^ \ \ 


V ^ — Stat. con. of Arutyunyan et. al. (1971) 


— Stat. con. of Arutyunyan et. al. (1971) A \ 

1.1 

jv —- Kep. seq. of Arutyunyan et. al. (1971) 

1.1 

— Kep. seq. of Arutyunyan et. al. (1971) A \' 


10^ 10^ 10^^ 10^^ 1.0 1.5 2.0 3.0 5.0 

Pc [ g/cm^] 3?eq [ 10^ km ] 


Fig. 3. Mass versus central density relation (left panel) and mass versus equatorial radius for 
general relativistic WDs using the Chandrasekhar EOS with p.= 2 for the static and the Keplerian 
sequence of this work and that of Arutyunyan et al. (1971). 


In Fig. [U the mass versus the central density (left panel) and the mass versus 
the equatorial radius (right panel) are shown for both Newtonian and GR WDs 
using Hartle’s formalism. As we can see in GR the maximum mass is reached at 
finite central density and its value is smaller than the Newtonian counterpart. Due 
to this fact in GR RWDs are subjected to the axisymmetric secular instabilities. 
Unlike GR WDs, Newtonian WDs reach their maximum mass at pc —>■ oo, hence 
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there is no axisymmetric instability in Newtonian WDs. The shaded region for 
GR WDs is the zone where uniformly rotating stable RWDs exist. One can also 
draw the stability region for the Newtonian WDs. Then this region will be limited 
from above by the mass shedding (the Keplerian) sequence and from below by 
the static sequence. Although in principle there is no limit from the right for the 
Chandrasekhar EOS because it depends only on the average composition and 
therefore no self-consistent inverse /3-decay instability line can be constructed for 
different chemical compositions, we assume here to guide the eye the critical inverse 
/3-decay density of "^He, 1.39 x 10^^ g cm“^. Heavier compositions have lower critical 
densities. 

Fig. [2] shows the eccentricity versus the central density (left panel) and the 
normalized quadrupole moment versus the central density (right panel). With the 
increasing central density the eccentricity and quadrupole moment start to decrease, 
which shows that the system becomes less oblate. As one can see owing to the 
general relativistic effects WDs in GR have smaller eccentricity and quadrupole 
moment than the Newtonian counterparts. 

Up to our knowledge, the only previous work on RWDs within GR is the one 
of Arutyunyan et. al.^^ A method to compute RWDs configurations accurate up to 
second order in the angular velocity of a star D was developed by two of the authors 
(see Ref. 22 for details), independently of the work of Hartle(1967.)^^ In Arutyunyan 
et. al.,^^ RWDs were computed for the Ghandrasekhar EOS with /i = 2. 

In Fig. [3] we show the mass-central density relation (left panel) and the mass- 
radius relation (right panel) obtained with their method with the ones constructed 
in this work for the same EOS. We note here that the results are different even at 
the level of static configurations, and since the methods are based on construction 
of rotating configurations from seed static ones, those differences extrapolate to 
the corresponding rotating objects. This fact is to be added to the possible addi¬ 
tional difference arising from the different way of approaching the order in the 
approximation scheme. The differences between the two equilibrium configurations 
are evident. 


Table 1. Maximum rotating mass of WDs in literature. 


Treatment/EOS 


References 

Newtonian/Chandrasekhar /i = 2 

1.474 

Anand (1965)^“ 

Newtonian/Polytrope n = 3 

1.487 

Roxburgh (1965)^® 

Post-Newtonian/Chandrasekhar /i = 2 

1.482 

Roxburgh and Durney (1966)^*^ 

GR/Chandrasekhar /i = 2 

1.478 

Arutyunyan et al. (1971)^^ 


Turning now to the problem of the maximum mass of a RWD, in Table [T] we 
present the previous results obtained in Newtonian, Post-Newtonian approach and 
GR by several authors. Depending on their method, approach, treatment, theory 
and numerical codes the authors showed different results. These maximum masses of 
RWDs are to be compared with the ones found in this work using the Ghandrasekhar 
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EOS /i = 2. In Newton’s gravity the maximum rotating mass is 1.534Mq and in 
GR is I. 5 I 6 M 0 . 

The minimum rotation period of WDs is obtained for a configuration rotating 
at the Keplerian angular velocity, at the critical inverse /3-decay density; i.e., this 
is the configuration lying at the crossing point between the mass-shedding and 
inverse /3-decay boundaries, see Fig. [T] We found 0.28 s for both Newtonian and 
general relativistic white dwarfs. It should be stressed that up to our knowledge the 
minimum rotation period has never been estimated for uniformly rotating white 
dwarfs. However Arutyunyan et al. (1971) showed their results for a range of central 
densities covering the range of interest of our analysis. Thus, we have interpolated 
their numerical values of the rotation period of WDs in the Keplerian sequence and 
calculated the precise value at the inverse /3-decay threshold for ^ = 2. We thus 
obtained minimum period ^0.31 s in agreement with our results. 

4. Conclusions 

We have constructed solutions of the Newtonian equilibrium equations for RWDs 
accurate up to order following the procedure of Hartle(1967).^^ In Figs. [I] [5] we 
have compared these Newtonian configurations with general relativistic RWDs for 
the Chandrasekhar EOS with p, = 2. We can see clearly the differences between the 
two mass-density relations toward the high density region, as expected. The most 
remarkable difference is the existence of axisymmetric instability boundary in the 
general relativistic case, absent in its Newtonian counterpart. 

We have shown, in particular, that the maximum mass of rotating WDs is stable 
against secular and dynamical instabilities and determined the minimum rotation 
period of rotating WDs both for the Newtonian and general relativistic RWDs. We 
have compared and contrasted the main differences between the Newtonian and 
general relativistic RWDs and we have shown that Newton’s gravity is not enough 
to properly investigate the stability and to calculate the maximum masses (rotating 
and non-rotating) of WDs. 

The results presented here open the way to a more general description of a 
rotating WD and are relevant both for the theory of delayed type la supernova 
explosions^^ as well as for the white dwarf model of Soft Gamma-Ray Repeaters 
and Anomalous X-Ray Pulsars. 
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